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The methods of non-equihbrium quantum field theory are used to investigate the possibifity of 
representing dissipation in the equation of motion for the expectation value of a scalar field by a 
friction term, such as is commonly included in phenomenological inflaton equations of motion. A 
sequence of approximations is exhibited which reduces the non-equilibrium theory to a set of local 
evolution equations. However, the adiabatic solution to these evolution equations which is needed to 
obtain a local equation of motion for the expectation value is not well defined; nor, therefore, is the 
friction coefficient. Thus, a non-equilibrium treatment is essential, even for a system that remains 
close to thermal equilibrium, and the formalism developed here provides one means of achieving this 
, numerically. 
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C/^ ; I. INTRODUCTION 

oo ; 

CNJ , Inflationary universe scenarios frequently refer to an inflaton equation of motion of the form 

+ 3i/</) + ?7((^)0 + ycff(<^) =0. (1.1) 

>: 

. The friction term ri{(j))(f> represents one mechanism through which the inflationary universe might be reheated, con- 
' verting the inflaton energy into particles ^ || , though the process of parametric amplification ^, ^ is in some 
. cases a more important effect. Moreover, it was pointed out long ago by Moss ^ that a frictional effect large enough 
' to maintain a thermal bath of particles during the inflationary era would have a significant effect both on the in- 
[ flaton evolution and on the primordial perturbation spectrum, alleviating the fine-tuning problems that afflict many 
inflationary models. More recently, this idea has been extensively investigated by Berera and co-workers |, p^ . 
While some frictional effect is to be expected in any interacting field theory, it is by no mea ns clear at a fundamental 
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. level that this effect can properly be represented by a local equation of motion such as ( |l.lD , and the purpose of the 

I ' present paper is to study this question in some detail. (A brief account of our findings has been given in Ref. |pT[.) 

Oh, To be fairly generic, we may consider a scalar field <&, which couples to other scalar fields and to fermions, denoted 

^ ' collectively by ^ and tp respectively. In Minkowski spacetime, a schematic Lagrangian density for this collection of 

'"^ ' fields is 



C = - im4$^ - iA$* - - .92*^^ + AC{C, V) (1-2) 



H I where AC is that part of the Lagrangian that does not depend on <I>. (The considerations developed here can probably 
be extended with little difficulty to a theory in which <I> also couples to gauge fields. They are essentially unchanged in 



a Robertson- Walker spacetime, as we discuss in section [X,) On splitting $(t, x) into its expectation value 0(t), which 



we take to be spatially homogeneous, and a fluctuation field ip{t, x) with {ip{t, x)) — 0, we find that the equation of 
motion for (^(<) is 

4> + m|0 -t- iA^3 ^ 1 A</,(^2^ + 310(^2^ + 52(^^) + lX{^3) = (1.3) 

while at tree level, the quantum fields ip, ^ and tp acquire effective masses 

m2(t) to| + iA02(i) m|(t) = TO^ +.gi(/)2(t) m^{t) ^ + g2(l>it) (1.4) 

which depend on time through (l){t). Dissipative effects in the equation of motion ( |1.3| ) can arise from the final 
term ^X{(p^), but we concentrate here on those arising from the other expectation values, which are quadratic in the 
quantum fields. Generically, the equation of motion is 

4>{t) + mUit) + g{x(t, xf)(t>(t) + ■ • • = (1.5) 
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where x may stand for Lp, the fluctuating part of $, for another scalar field ^ or for a fermionic field -0, and has a 
0(i)-dependent mass of the form (1.4). In general, there will be a sum of such contributions from the fields appearing 
in a specific model, but our considerations are substantially model-independent. 

Superficially, it is quite plausible that a system which remains fairly close to equilibrium can be treated by using 
equilibrium statistical mechanics. In section ^ we review the relevant dissipative mechanisms and the results for 
the friction coefficient ?7((/') which have been obtained on this basis; we also indicate why this approach might be 
questioned, when applied to a system that is not maintained in exact thermal equilibrium. Sections III-^ develop 
a strategy for tackling the problem in non-equilibrium scalar field theory, exhibiting the sequence of approximations 
necessary to obtain a local equation of motion of the form (1.1). We find that this local equation of motion contains 
coefficients whose values must be found from the solution of auxiliary kinetic equations. The resulting set of local 
evolution equations is suitable for a numerical solution. However, we investigate in section VI whether an adiabatic 
treatment of these evolution equations leads to a well defined friction coefficient rj{(t)). We find that it does not, 
and section VI] exhibits numerical evidence that non-equilibrium effects may be quantitatively quite significant. 



Section VIII briefly discusses frictional effects due to fermions, showing that the formal situation is quite similar to 
that developed in detail for scalar flelds. Finally, our principal conclusions are summarized in section where we 
also discuss the extent to which they depend on the methods of approximation we have utilized. 



II. MECHANISMS OF ENERGY TRANSFER AND DISSIPATION 



Roughly speaking, we can identify two mechanisms through which energy may be transferred between the classical 
field <f> and the system of x particles. One is the creation of new particles, which we will refer to as type-I. The other, 
to which we refer as type-II, involves changing the energies of particles that are already present. In the following 
subsections, we review the arguments which purport to derive from these two mechanisms friction terms in the 
equation of motion for 0. 



A. Particle creation 



A crude argument given by Morikawa and Sasaki takes x iii the first instance to be a free field, except 
that the i nteraction with leads to a time-dependent effective mass. Thus, the energy of a single x particle is 
ujk{t) = y/k"^ + TTi^ + gfj)^ (t) ■ According to standard arguments, this field can be represented as 



(2^ 



e'^--fk{t-,t)ak{t)+e-^'^-'^m-t)al{t) 



where the mode functions fk{t;t) are solutions of the equation 
They satisfy the Wronskian condition 



(2.1) 



(2.2) 



(2.3) 



where the overdot indicates differentiation with respect to the first argument, t. Of the many complex functions 
that obey these two equations, we select a family, parametrized by the reference time t, by imposing the boundary 
conditions 



These mode functions can be expressed as 



-1/2 



fk{i;i) = -i [LJk{i)/2] 



1/2 



fkit;i) = [2nkit;i)] '^'exp 
where the frequency r2fc(t;i) is the solution of 



-i J Vtk{t']i)dt' 



IVlk 3 ^2? 2 _ , ,2 



(2.4) 



(2.5) 



(2.6) 
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subject to the boundary conditions flkit;t) — ujk{t) and Ctk{t;t) — (lk{t;t) = 0. At times close to t, therefore, fk{t;t) 
is the positive-energy solution, 



fk{t;t) « [2u;k{i)] 



'1/2 



exp 



-iujk{i){t - i)] 



(2.7) 



but at other times the frequency flk{t;t) does not necessarily correspond to a single-particle energy. Mode functions 
referred to different reference times are related by a Bogoliubov transformation. To be concrete, let /°(i) = fk{t;0). 
Then fk{t;t) is a Hnear combination of fk{t) and fk*{t), say fk{t;t) = A{i)f^{t) + B{t)f^*{t). It is a simple exercise 
using ( p^ ) and (2.4) to find the BogoHubov coefficients A{i) and B{t) and thus the dependence on i of fk{t]i), which 
is 



I UJk{t) 



(2.8) 



The creation and annihilation operators in have the commutation relation [afc(£), a|,, (£)] = (27r)'^(5(A; — fe'). 

Their dependence on the reference time t follows from the fact that x{^^ itself is independent of t; using the fact 
that fk{t; t) depends only on fc = |fc|, we find 



diak{i) = -iujk{i)ak{i) 



luJkii) t a\ 

^a_j,(r). 

2a;fe(i) 



(2.9) 



With this formalism in hand, we evaluate the expectation value in (L5) by choosing the reference time i to be t. 
For a spatially homogeneous state, the result is 



{X{t,xf)) 



(Pk 



{2^f2ujk{t) 

where the functions nk{t) and Vkit) are defined by 



[l + 2nk{t) + 2Reuk{t)], 



{al{t)ak'{t)) - {2TTf5{k-k')nk{t) 
{ak{t)ak'{t)) ^ {2T:f5{k + k')uk{t). 



(2.10) 



(2.11) 
(2.12) 



The functions Ukit) can be interpreted as occupation numbers for the single-particle modes fk{t]t) when t is close to 
t, while i^kit) measures the off-diagonality of the density matrix in the representation defined by these modes. The 
time dependence of these functions is easily established from (2^), and we find 



dtukit) = '^Rei.kit) 

OJk(t) 

dtVk{i) = -2iujk(t)vk{t) + 



2L0k{t) 



[l + 2nk{t)] 



(2.13) 
(2.14) 



These evolution equations are valid only in the approximation that x(t,a;) is a free field. It is argued in that 
interactions will give rise to an imaginary part of the x self-energy, which can be taken into account by replacing 
LOk{t) with ujk{t) — i^k{t)- With the further assumption that ujk{t), ^k{i) and nk{t) change negligibly on a time scale 
of order Tkit)^"^, integration of (2.14) yields 



Rez/fc(t) 



UJk{t)Tk{t)uJk{t) 



[l + 2nk{t)] 



if Ffe is negligible by comparison with Uk- Since ujk 
equation of motion (1.5), with 



2[u;l{t)+Tl{t)] 

LOk, we identify a contribution of the form rji 



Tk{t) 



{2nrcukmLolit)+rlitW 



[l + 2nk(t)]. 



(2.15) 
in the 

(2.16) 
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B. Changing particle energies 



To the extent that the x particles can be considered as free, the energy of this system of particles can be estimated 

as 

(2^ 



-^W - I 7KZ^z^k{t)nk{t). (2.17) 



Clearly, this energy is altered by a change in the classical field 4){t), which causes a change in the single-particle 
energies ujk{t), and this provides a mechanism whereby energy may be exchanged. However, a frictional term in 
the equation of motion for arises only from an irreversible transfer of energy. This is a secondary effect, brought 
about by the fact that a change in ujk{t) alters the scattering and decay rates of x particles, and hence affects the 
evolution of the occupation numbers nk{t). A crude argument for estimating the friction coefhcient that arises from 
this mechanism was given by Hosoya and Sakagami [ p^ , who take the time evolution of nk{t) to be governed by a 
phenomenological kinetic equation of the form 

dtukit) = ~2Tk [nk{t) - nl^{t)] , (2.18) 

which is sometimes referred to as the relaxation-time approximation to the Boltzmann equation. It is implied that 
the system of particles is always close to an equilibrium state characterized by a fixed temperature f3~^, with the 
distribution function 

<^(t) = l/[exp(/3c^fc(i))-l]. (2.19) 



If we once more assume that i\^{t) changes slowly on a time scale ^, then (2.18) can be approximately integrated 
to yield nk{t) = rij^{t) + 5nk{t), with 

Snk{t) « -^{2Tk)-'dtn''i{t) 

« /3(2rfc)-i^^fe(t)n^'i(i)[l + n°^(t)] (2.20) 

and we obtain a contribution to the equation of motion ( |1.5[ ) of the form 7711(0)'^, with 

(3g^^t) f d^k n^^it) [l + n°q(0] 

C. Linear response theory 

The crude arguments given above can in some respects be improved by restricting attention to a situation in which 
the system of x particles is in thermal equilibrium, apart from a small time-dependent perturbation which is treated 
to linear order. Different versions of this treatment have been given by Hosoya and Sakagami [ p^ , by Morikawa and 
Sasaki jlj, [ij] and more recently by Berera, Gleiser and Ramos [|l5[ |l^, |l^. It is necessary to suppose that over a 
sufficiently long period of time, the classical field (t>{t) can be decomposed as (/)(t) ~ (t)Q + S(t){t), where (f>o is constant 
and 64)(t) is small. Correspondingly, the Hamiltonian for x will be decomposed as 

Hix,t) = Hoix)+g<f'oSHt) j d'xx\t,x) + 0{5<j>^), (2.22) 

the second term being treated as a perturbation. To linear order in this perturbation, the standard Kubo formalism 
then yields 

(X^(t,a;))«(x2(t,a;)r + z#0 f dt'5cf>{t') f d^x' {[x^if ,x'),x^{t,x)]r (2.23) 

where ""^ denotes the thermal average in the equilibrium ensemble determined by Hq{x)- Although -ffo(x) has no 
explicit time dependence, it still contains interactions, and the remaining thermal averages can be computed sys- 
tematically in perturbation theory. To obtain a well-defined answer for the linear response term in ( 2.23| ), it proves 



necessary to effect a partial resummation of the x propagator, in particular identifying the thermal width as the 
imaginary part of a suitable self-energy (as discussed in more detail below). As first pointed out in Ref. jlj], the net 
result is a contribution to t he eq uation of motion (|l.5|) of the form rj{(j))4>, with ?7((/)) « rji{(f)Q) + r]u{(t>Q), provided that 



Tk ^ ujk and that nk(t) in ( 2.16 ) and n°'i{t) in ( p.2lD are identified as the constant equilibrium distribution associated 
with iJo(x)- As in the previous calculations, it is also necessary to assume that 6(j){t') varies slowly on a time scale of 
order F^T^ and can be approximated as S(j){t') « (j){t){t' — t). 
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D. Open questions 



Each of the calculations outlined above has its own deficiencies. The type-I analysis of section [I A is valid arbitrarily 



far from equilibrium, but treats interactions in an incomplete and ad hoc manner. The type-II analysis of section 



[IB is explicitly restricted to states very close to local eq uilibr ium and rests on a kinetic equation which is little 



more tha n a g uess. The linear response treatment of section II C is much more systematic, insofar as the expectation 



value in (2.23) can in principle be evaluated at any desired order of equilibrium perturbation theory. However, this 
calculation depends in an essential way on analytic properties of thermal Green functions that are meaningful only 
in a state of exact thermal equilibrium. Specifically, the Wightman functions [t — t';x — x') = {x{t, x)x{t' , x'))'^'^ 
and G^(i — t';x ~ x') = {x{t' , x')x{t, x))°'^ have Fourier transforms which satisfy the KMS condition G^{ko,k) = 
exp(/3A;o)G^ (fco, fc) (see, for example, Rcf. fl^). In a state which departs even slightly from thermal equilibrium, 
neither the temporal Fourier transform nor the KMS condition has any meaning. 

The primary question addressed in the remainder of this paper is whether the apparently plausible result of linear 
response theory is recovered, for a non-equilibrium system, in the limit of slow time evolution. This will be possible, 
at best, only if we have a local approximation to expectation values such as (x^(i, x)), which are inherently nonlocal 



in time (as evidenced by ( p.2S ) even in the linear response approximation) , and devising such an approximation is 



the key feature of the analysis that follows. 

An important issue that is not addressed in this paper arises from the expression (2.21). This contribution to 
the friction coefficient depends inversely on the thermal width F^ , which is typically of the order of the square of a 
coupling constant, and appears to call into question the reliability of perturbation theory as applied to this problem. 
In fact, this is typical of expressions obtained in the application of linear response theory to the estimation of transport 
coefficients, and it is known that infinite classes of diagrams contribute at each order of perturbation theory |jl^, |2^, 
How the requisite resummation might be effected for a nonequilibrium system is beyond the scope of this paper. 

III. LOCAL APPROXIMATION FOR SELF-ENERGIES 

The exact two-point functions for the quantum field x(i, x) will be denoted in a standard notation by 

Gab{t,t';k) = J d^xe'''-''Gabit,x;t\0), (3.1) 

(for a,b — 1,2) where 

x,t,x)^^ ^^^^^ ^^j^^^^ J . [6.Z) 

Self-energies Yiab{t,t']k) can be defined by the Dyson-Schwinger equations 

Gab{t, t'-k)^ g^J> {t, t';k)-i f dt"dt"'gZ> {t, t"; fc)S,rf(t", t'"; k)Gdb{t"' , t'; fe), (3.3) 



( F ) 

in which the frcc-ficld propagators g^b {'^^ t'] k) are solutions of the equations 

Pi^)(t, duk)g^\t, t'- k) = 5(F) (i, fc)p(P(t', dr-k) = ^i5ab5{t - t'), (3.4) 
where the differential operator P'^^^(i, dt; k) is given by 

_(d^ + e + m\t) \ 

\ -d'^ -e -m'it)) ^•^■^> 

with m?{t) — m? -\- g(f>'^(t). The form of the self-energy matrix Y.ab is constrained by some general considerations. 



First, the full propagators Gabit, t'; k) defined by ( |3.l[ ) and (3^) have the properties 



Gab{t,t'-k)=Gba{t',t;k) (3.6) 

Gti(t,i';fc) =G22(t,<';fe), Gl^{t,t'-k)^G2i{t,t'-k). (3.7) 
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(a) (b) (c) (d) 

FIG. 1: Diagrams contributing to the 2-point function of X^"^ theory. Diagram (a) gives a local contribution to the self-energy; 
diagram (c) gives rise to the principal dissipative effects; diagram (d) represents the contribution of Mab- 



Second, causality requires the integrand in ( ^^ ) to vanish if either t" or t'" is larger than both t and t' . From these 
observations, it is not hard to show that the self-energy matrix has the general form 



E{t,t';k) = j:^{t;k)S{t-t') + j:>{t,t';k)e{t-t') + j:<{t,t';k)e{t' ~t) (3.8) 



with 



T.<{t,t';k) = (^j^*[/J^ -AU/J)) ■ ^^-^^^ 

In the local part S'^, the functions ak{t) and (3k{t) are real. If x has, for example, a quartic self-coupling x^, or a 
biquadratic coupling y^^^ to another scalar field ^, then these couplings generate Feynman diagrams (for example, 
diagram (a) in figure |^) that are manifestly local; such diagrams contribute to /3fc (t) , but not to ak (t) . More generally, 
it will be advantageous to extract local contributions from non-loc al d iagrams, and it is easily seen that a contribution 



to Ak{t,t') of the form [Pk{t) — iafc(i)]^(t — t') is equivalent to (3^), provided that we identify 6'(0) = \. We take 
Q!fe(i) = for now, but will later use afc(0 to denote a local contribution to the imaginary part of Afc(t, t'). 

In fact, the goal of this section is to develop a plausible ansatz for approximating the whole self-energy matrix as 
a local quantity, containing terms proportional only to S{t — t') and dt5{t — t'). To this end, it is convenient to deal 
with the commutator and anticommutator (or correlation) functions 



p{t,t'-k) = I J d^xe^''-^{[x{t,x),xit'M) 

= i[G2i{t,t'-k)-Gi2{t,t'-k)] (3.12) 
C{t,t'-k) = i Jd^xe^'^-^{{x{t,x),x{t\0)}) 

= i[G2i(t,t';fe) + Gi2(t,i';fc)] (3.13) 
from which all the Gab can be constructed. For these functions, the Dyson-Schwinger equations (|3.3[) together with 



(3.4) imply the equations of motion (see, for example, Ref. |2S 



[d^ +ujl{t)]p{t,t';k) = - / dt"Ysp{t,t"-k)p{t\t'-k) (3.14) 



(3.15) 



[d^^ +u}l{t)]G{t,t']k) = - / dt"T.p{t,t";k)C{t",t';k)+ f dt" Y.cit,t";k)p{t" ,t;k), 

Jo Jo 



with 



cuKt) = e+m'^{t)+Pk{t) (3.16) 
Sp(t, t'; k) = [Ak{t, t') + Al{t, t')] e{t - t') - [Ak{t', t) + Al{t', t)] e{t' - t) (3.17) 
Ec(t, t'; k) = -li [Ak{t, t') - Al{t, t')] 0{t - t') - li [Au{t', t) - Al{t' , t)] d{t' - t). 

(3.18) 



7 



In a state of thermal equilibrium, for which p{t,t'; k) = p^'^{t — t'; k), one can define the spectral density 

p°^(a;,fc) = -i / dte^^'p'^^mk). (3.19) 



A strategy frequently adopted in the context of linear response theory is to assume that this spectral density can be 
approximated by the Breit-Wigner form 



p^\uj,k) 



(3.20) 



where ^Ik and are a quasiparticle energy and width to be abstracted from the real and imaginary parts of the 
self-energies. The numerical calculations of Ref. for a one-dimensional scalar field theory suggest that such an ap- 
proximation is reasonable (though certainly not exact) for a non-equilibrium state also. To realize this approximation 
in the present context, observe that the real-time commutator function is given by the Breit-Wigner approximation 
as 



p^^it-t'-k) 



^-ink(t-t') _ ginfc(t-t') 



-Tk\t-t'\ 



(3.21) 



For the non-equilibrium state, we adopt an ansatz which is a natural generalization of this expression, namely 
p{t,t';k) « p^^\t,t']k), where 



Qk{t,t') 



exp 



dt"Tk{t") 



0{t - 1') 



(3.22) 
(3.23) 



and fk {t; t) are the mode functions introduced in section II A , except that the quasiparticle energies oJk (t) may include 
loop corrections. It is easily checked that p^^\t, t'; k) is real and antisymmetric in its time arguments, as it should be. 
It also has the property dtp'^^\t,t'] k)\^,_^ = 1 as required by the canonical commutator. Finally, as can be verified 

from (2.8), it is independent of the reference time provided that Vk{t) is. The approximate version of the equation 
of motion (3.14) obeyed by p'^^\t,t'] k) is conveniently expressed by writing 



pW(t, t'- k) = p(°)>(t, t'- k)e{t - 1') + p'-°^<{t, t'- k)e{t' - 1), 

with p^°^<{t,t';k) = -p'^">{t',t;k). By differentiating ( |3^ , we find 

+ 2Tk{t)dt+ojl{t)+Tl{t)+rk{t)]p^''>{t,t';k) - 
'df - 2rkm +Lolit) +TUt) ~tk{t)] p^''^<{t,t';k) = 0, 
together with the equal-time conditions 



V^<(t,t';fe) 



(3.24) 

(3.25a) 
(3.25b) 

(3.26) 



Comparing these equations with (3.14), we see that they imply a corresponding ansatz for the self-energy ( |3.17D , or 
for the real part of Ak{t,t'). It is 



Re (t,i') = [Vk{t)Tk{t')f^ dt5{t - t' - e), 



(3.27) 



wh ere th e positive infinitesimal e is included to ensure that the delta function is satisfied inside the range of integration 
in ( ini ). 

Given a nonzero quasiparticle width Tk{t), the non-local self-energies may be expected to decay roughly exponen- 
tially with \t — t'\. For T,c{t, t'; k), which is symmetric in its time arguments, or, equivalently, for the imaginary part 
of Ak{t,t'), a suitable local ansatz is 



luiAf{t,t')^^ak{t)5{t~t'). 



With this approximation, the equation of motion (3.15) for the correlation function becomes 

d'i + 2Tk{t)dt + Lol{t) + Tl{t) + f fe(t)l {t, t'; k) - -ak{t)p^°\t, t'; k)9{t' - t). 



(3.28) 



(3.29) 
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Alternatively, defining 

C(°) (t, i'; k) = C'-">{t, t'- k)e{t - t') + C(°)< (t, t'- k)0{t' - t), 
with C(°)<(t,t';fc) = C(°)>(t',t; fc), this may be written as 

+ 2Tk{t)dt+ujl{t) +Tl{t) + tkit)\c^''>{t,t'-k) = 
'd^t+2Tu{t)dt+ujl{t)+Tl{t) + tk{t)\ C^°^<{t,t';k) = -ak{t)p^°^<{t,t';k) 
together with the equal-time condition 

dtC^°>{t,t';k) ^ dtC^°^<{t,t';k) 



(3.30) 

(3.31a) 
(3.31b) 

(3.32) 



Finally, we assemble the approximate correlation and commutator functions into a single complex function, 

C7(")>(t,t',fc) ~ iip(")>(t,t';fc) = hkit,t') (3.33a) 
C7(°)<(t,t',fe) + iip(°)<(t,t';fe) = hk{t',t). (3.33b) 



Then, denoting by gab{i, t'; k) the approximation to the propagator matrix (3.1) that embodies the ansdtze ( p. 27 ) and 
(3.28), this matrix can be summarized by 



7a6(i, t'- k) = hb{t, t'- k)0{t - t') + ha{t', t- k)e{t' - t), 



(3.34) 



with hi{t,t';k) — hk{t,t') and h2{t,t';k) — hl{t,t'). The approximate equations of motion ( 3.25| ) and (3.31) are 

+ 2r kirn + Lol{t) + rlit) + r kit)] hk{t,t') = o (3.35) 



[d?+LuUt) + Tlit) - iakit)] hk{t',t) + hVkim + rfc(i) + iak{t)\ Kit' ,t) = 



while the equal-time conditions ( |3.26 ) and ( 3.32 ) become 



dt[hk{t,t')-hk{t'Mt'=t^~^ dt[hl{t,t')~hk{t'Mt'=t^^- 



(3.36) 



(3.37) 



At this point, the ansdtze ( ^.27 ) and (3.28) have yielded local equations of motion for the approximate two-point 
functions, but the functions ak{t) and Tk{t) that appear in these equations and the function I3k{t) that appears in 
the single-particle energy (3.16) are unknown. Given a specific Lagrangian, approximations to the self-energies Sofc 
can be obtained — for example, from some version of perturbation theory. The function I3k{t) can be identified from 
the local part of Eafc, bu t a p rescri ption is needed for extracting from the calculated Ep and Sc local contributions 
of the kind indicated in (3.27) and (3.28) so as to identify Vk(t) and ak(t). This issue will be addressed in section 0. 



IV. LOCAL KINETIC EQUATIONS 



The approximation summarized by (3.34) - (3.37) is equivalent to that obtained from somewhat different consider- 
ations in Ref . ||2^ . There it was shown that the general solution to the local equations of motion can be expressed in 
the form 



hk{t,t') ^ iexp 



dt"Tk{t") 



hkit,t') 



hit, t') = [1 + Qk{t'; t)] h{t; i)f*{t': t) + [-1 + Q*(t'; t)] /,*(t; i)h{t'- 1). 
The function Qk{t;t) is a solution of the equation 



dt + 2Tk{t) + 2 



mi) 



[dt + 2Vk{t)]Qk{t]i)^2iak{t) 



(4.1) 
(4.2) 

(4.3) 
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subject to the constraint 



{dt + 2Tk){Qk + Ql) + ^{hft)-\Qk - Ql) = 0, 



(4.4) 



which is preserved by (4.2). For our present purpose, it is useful to observe that Qk{t]i) can be decomposed as 

gfc(t;£)^Q«(t;f) + 4MQf 



where Q')^\t; i) is real. If Q)^' and Q)^' are taken to obey the first-order equations 

[dt + 2Tk{t)] Q[^\t; i) = 2fk{t; 1)^(1; i)ak{t) 



l(2) 



[dt + 2Tk{t)]Q'^^\t;t) 



(t;£)afe(t), 



(4.5) 



(4.6) 
(4.7) 



then both (O) and (44) are satisfied and the decomposition (4^) is unique. 
The propagator function /ifc(t, t') is now given by 

hkM = [l + Q«(t';t)l h{t;i)rk{t';t)+\-l + Q^^\t'-i)\ r,{t;i)fk{t' -J) 



+ Qi'\t'; i)fk{t; i)h{t'; i) + Q)^^*{t'; i)f^{t; i)f^{t': t) 



(4.8) 



It must be independent of the reference time and this determines the dependence on i of the functions Q^^\t]i). 
We find 



In the equation of motion (1.5), our approximation to {x (t,x)) is 



(2^ 



hk{t,t) 



and this is conveniently evaluated by choosing the reference time t to be the time t of interest: 

hk{t,t) = Q[}\t;t) + ReQ^^\t;t). 

In fact, let us define 



nk{t) 



We have 



hk it,t) = [2ujk (<)] Ml + 2nfc [t) + 2Re Vk [t)] 



(4.9) 
(4.10) 

(4.11) 

(4.12) 

(4.13) 
(4.14) 

(4.15) 



and will loosely identify the functions nk{t) and Vkit) with the quantities denoted by the same symbols in section ||. 
By combining (4.6) and ( [1.7] ) with (4.9) and ( 4.1Cl| ), we obtain the evolution equations 



dtnkd) = ^^^Tk{t)[l + 2nk{t)] + '^Revk{t) 



2ujk{t) 



Uk{t) 



dtuk{t) - ~2i[LOk{t)-iTk{t)]vk{t)-^^,+^^,[l + 2nk{t)]. 



(4.16) 
(4.17) 



2uok{t) 2uJk{t) 

These are clearly generalizations of the free- field equations ( 2.13| ) and ( 2.14 ), the extra terms involving the functions 
ak{t) and Tk{t) whose exact meanings are explored further below. On the other hand, the tentative identification 

ak{t) ~ 2iUkit)Tkit) [1 + 2n^^(t)] (4.18) 



(which we will later find to be too naive) brings (4.16) into the form 

dtUkit) ^ -2Tkit) [nk{t) ~ nl^it)] 



LUk{t) 



Revk{t), 



(4.19) 



which is a generalization of the kinetic equation (2.18) including a source term to account for particle creation. 
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DETERMINATION OF LOCAL SELF-ENERGIES 



To give substance to the kinetic equations ( 4.16| ) and (4.17), we need a concrete method of determining the functions 
f^kit) (which appears in the quasiparticle energy ( |3.16| )), 0.^(1) and Tf;{t). A prescription for doing this was given in 
p3| ; here we describ e a re finement of th at pr escript ion w hich is convenient for the problem at hand. The approximate 
two-point functions ( 3.34 ) which solve ( 3.35 ) and ( 3.36| ) are the exact propagators of an approximate theory defined 
by the closed-time-path action 



'5'cTp(Xi,X2) = -5 J dt J d^xXa{t,x)Vab{t,dt,V)xb{t,x), 



where the differential operator V is given, after a spatial Fourier transformation, by 

V = 



-2rkit)dt-tk.{t) + iak{t) 



2Tkit)dt+Tkit)+iakit) 



(5.1) 



(5.2) 



(An effective action having essentially this structure also describes an open system, coupled to environmental degrees 
of freedom, which can be integrated out by the Feynman- Vernon influence functional method ||2j, 25, Here, one 
may think of a single field mode having an environment that consists of all the other modes, but this environment 
is treated in a self-consistent manner, rather than being integrated out.) If the complete theory has the action S{x), 
and the corresponding closed-time-path action 5'ctp(xi,X2) = 'S'(xi) — S'(x2), then a partly resummed perturbation 
expansion can be defined by writing 



(5.3) 



<5'CTP(X1,X2) = 'S'cTp(Xl,X2) + AS'cTp(Xl,X2) 

and treating AS'ctp as the perturbation. Included in A^ctp is the counterterm i J dt cfi x Xa-MatXb, with 

M{t,dt;k) = V{t,duk)~V^''\t,duk) 

fik{t) + Tl[t) - iak{t) 2Tk{t)dt + tk{t) + iak{t) 
-2rk{t)dt - tk{t) + iak{t) -I3k{t) - Tl{t) - iak{t) 

which accounts for the difference between gab and .g^f, . In the context of this expansion, self-energies are defined by 
replacing g^^^ in the Dyson- Schwinger equation (3^) with gab- They have the form 



(5.4) 



Y.ab{t, t'; k) - Y}°°^{t, t'- k) - Mab{t, dt;k)S{t - t'), 



(5.5) 



where y}°°p consists of loop diagrams in which the propagators are gab- 

We would like to choose afc(i), (3k{t) and Tk(t) in such a way that the propagators gab{t,t' ; k) approximate the 
exact two-point functions as closely as possible. Loosely speaking, this means making the self-energies (5.5) as small 



as possible. More precisely, it is necessary to obtain a local approximation to , which can be cancelled by an 
appropriate choice of TWah- To this end, suppose that t and t' are both close to the reference time t. We define 
time-translation-invariant propagators which approximate gab{t,t'; k) in this region by introducing the function 



hk{i,T) 



.-PfcT 



2a) A: 



{[1 + nk{i) + i^kme-'""-^ + K(£) + i^liW'^} , 



(5.6) 



where t = t — t' , ujk — ujk{t) and — ^k{t\ This function is obtained from the one defined in (4.1) by taking 
hk{t,T) = hk{t + T^t), using the approximation (2.7) for fk{t;t), and replacing Tk{t") with Tk in the prefactor. The 
approximate propagators are then given by 



gab{i, t; k) ^ hb{i, t; fc)6'(T) + ha{i, -t; k)0{-T), 



(5.7) 



where, as in ( 3.34 ), we use the notation hi{t,T;k) = hk(t,T) and h2{t,T;k) — h*f,[t,T). 

By replacing gab with gab in the diagrams that constitute y}"^"^ , we obtain a time-translation invariant approximation 
to these self-energies, r; fc), valid when t and t' are both close to i, which can be used to determine the functions 

ak{i), fik{i) and Tk{i). After a Fourier transform on r, our approximation to the right-hand side of (^.E|) is 



duj 
2^ 



t}°°''(i-u:,k) - Mab{t,-iu^-k) 



(5.8) 
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Although E'°°P(i,T;fc) is time-translation invariant, it is not, in general, a distribution concentrated at r = 0. In 
Fourier transfo rmed language, S'°°p(£; cj, fc) is a non-linear function of l o, w hile Mabit, —ii^', k) (which results from 
replacing dt in ( ^.4[ ) with —iu) is linear in uj. Therefore, the integrand in ( |5.8| ) cannot be made to vanish for all values 
of Lo. A reasonable prescription for determining ak{t), I3k{t) and Tk{t) is to demand that 



(5.9) 



To the extent that the considerations of section III are valid, this ensures that gab{t, t'; k) are the propagators for free 
quasiparticles whose energies and widths are approximat ely t hose determined by the peaks of the true non-equilibrium 
spectral density. Explicitly, the prescription implied by (5.9) is 



ak{i) 
Tkii) 
Pk{i) 



t'°°^{i-Qu.k) + t'°°^k -^k,k) 
^^21 ^{i; <^k,k)- f^2i^{i; -LOk, k) 



I 

ReSi°°P(f;^fc,fc)-r2(£). 



(5.10) 
(5.11) 
(5.12) 



VI. FRICTION TERMS IN THE EQUATION OF MOTION 



The nature and purpose of the approximations we have introduced so far can usefully be summarized as follows. 
By partitioning the closed-time-path action as in (b.3| ), we obtain a reorganized perturbation theory in which the 
unperturbed propagators are the gab defined by (|3.34|) - (3.37). Summed to all orders, this perturbation theory would 



(presumably) be equivalent to the usual expansion based on the free-particle propagators g^b , ^-^d this assertion 
is essentially independent of how we choose the functions Q!fe(i), (3k{t) and Tk(t) that enter the definition of gab- 
However, the reorganized perturbation theory cannot in practice be summed to all orders. Its utility rests on the 
possibility of making gab a better approximation than g^j, to the full two-point functions; in particular, we wish to 
estimate the expectation value {^^{t^x)) by retaining only the lowest-order term as indicated in (4.11). To do this, 
we choose ak{t), f3k{t) and Tk{t) to be local contributions to the true self-energies, which are thereby resummed in 
the reorganized perturbation expansion. For the purpose of extracting these local contributions (and only for this 
purpose) we introduced t he ti me- trans lation-invariant propagators gab m (5/7), which enabled us to formulate the 
prescription recorded in (5.1C) - ( 5.12| ). Of course, this prescription can be implemented only approximately, by 
evaluating 2'°°? to some finite order of perturbation theory. 



In principle, we are now in a position to evaluate the right-hand sides of the evoluti on en uations (4.16) and (4.17), to 
solve these equations for nk{t) and I'kit) and hence to estimate the expectation value (4.11) in which we are principally 



interested. There is, however, a practical difficulty. It is that the self-energy on the right-hand side of (5.11) is itself 
a function of F^, and this equation cannot be solved analytically to obtain a concrete expression for F^. A numerical 
solution is feasible, and this is no doubt the best way of estimating the time evolution, given a specific model. Here, 
though, we wish to investigate the circumstances under which dissipation might be represented by the frictional terms 
in the equation of motion exhibited in section H. To that end, we now introduce two further approximations. First, we 
take the limit Ffe ^ in the propagators (jab used to calculate Eob in ( 5.10| ) - (5.12). This is reasonably well justified 
in a weakly coupled theory, where F^ is of order g^, say, and Eat itself contains an ov eral l factor of g^. Second, we 
will also set Vk — in gab- A justification for this step will appear below [see equation (6.7)]. 

With these approximations, the propagators gab assume the form familiar from the equilibrium theory (see, for 
example, Ref. pSl). In particular, the temporal Fourier transform of 512 is 



g2iii;uj,k) 



i^k 



{[1 + nk{i)] S{uj -Uk) + nk{i)5{uj + Cok)] 



(6.1) 



As a standard exam ple, we consider in what follows the A<f>^ theory, whose Lagrangian density consists of the first 
three terms of (1.2), identifying x as the fiuctuation field and the coupling g as g = A/2. The 2-point functions 
for X contain, amongst others, the diagrams shown in Fig. |^. At 2-loop order, with ^21 given by (6.1), the only 
contribution to E21 comes from diagram (c). (The 1-loop diagram (6) contains products of 6 functions which cannot 
be simultaneously satisfied. It would give an off-shell contribution if we were to retain a non-zero width Ffe in ,921)- 
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Evaluating this diagram, we find from (5.10) and (5.11) 



ak{i) 



A2 



32(27r) 



(fki(fk2d^k 



,3 S{uki + '^k2 - '^fca - (^k)5{ki +k2-k3-k) 



X [(1 + rife J (1 + nfejnfc., + nk^nk^il + ^^fe;,)] 



A2 



64(27r) 



S{LOkt + LUk2 - l^ka - '^k)S{ki + fc2 - ^3 - fc) 



^k^ki'^k2^k3 



X [(1 + nfej(l + nk2)nk3 - fik^fik^il + ^^fe3 
where n^. — nki{t). In particular, the quantity 



(6.2) 
(6.3) 



Sk{i) 



Tkii) [l + 27ife( 



akji) 

2Lbk 

32(27r)''5 J LOkUJkiUJk^uJks 
X [nkifik^il + fik3){l + fik) - (1 + '^fci)(l + nk2)nk3fik] 



d^k^d'k^d^ks^^ 



^ks 



a)fe)(5(fci + ^2 - fcs - fe) 



(6.4) 



which appears in the kinetic equation (4.16) is precisely the 2-particle elastic scattering integral that ought to appear 
in a genuine Boltzmann equation. (More generally, since the A^'"' theory has no conserved particle number, inelastic 
contributions should also be expected, and these will indeed appear if we extend the evaluation of the self-energies 
beyond 2-loop order.) 



Wi th th ese explicit expressions for akit) and Tk(t), we now have concrete forms for the evolution equations (4.16) 



and (4.17). These evolution equations are local in time, but their solutions will be non-local. To extract contributions 
to {x^itjxY) which are proportional to 0(t), we must resort to some adiabatic approximation of the kind considered 
in section ||. To do this systematically, we rewrite the evolution equations as 



edtUkit) = dkit) -rk{t)[l + 2nk{t)]+e 



i^k{t) 

UJk{t) 



edtVk{t) 



-2i[ujk{t) - iTk{t)]vk{t) - c^k{t) 



Rcvkit) 

^ (^k{t) 
'Su^kit) 



[l + 2nfe(t)] 



(6.5) 
(6.6) 



where ak{t) — ak{t)/2L0k{t) and we have introduced a formal expansion parameter e multiplying terms with time 
derivatives. By expanding in powers of e and finally setting e — 1, we generate expansio ns of nk{t) and Vkit) in 
powers of the time derivatives of (j){t). On substituting these expansions in the expression ( 4.1l| ) for {x^(t,x)), the 
next-to- leading terms, proportional to i^, yield an estimate of the friction coefficient ri{(j)). 

At leading order, ( |6.5D reduces to the equation Sk{t) — whose solution is well-known to be the Bose-Einstein 
distribution nk{t) = nj^{t) for some temperature (We do not allow for a non-zero chemical potential, because 

the inelastic contributions to Sk(t) expected at higher orders of p erturbation theory would constrain the equilibrium 
chemical potential to vanish.) The corresponding solution of (|6.6[) is 



where a^'^ and T'^^ are obtained from 



l.°^Ht)^^at'/2[uJk{t)-^Tl\t)] 
and (|6.3|) by setting Uk — i 



smaller than by a factor of order Tk/uJk, so setting Vk 
approximation. 



(6.7) 

We see that, near equilibrium, Vk is 
in the calculation of the self-energies should be a fair 



At next-to-leading order, we set Uk 



eSuk and Vk 



e6i>k- The linear-response approximation to ??((/)) 



the sum of (2.16) and (2.21)] can be recovered only at the expense of further approximations. The first is to replace 



Tkin) with r^"^ = rfe(n^) and similarly for afc(n): that is, to ignore the fluctuations in these quantities brought 
about by fluctuations in nk{t). The second (when calculating dty"^^) is to take the time dependence of r|,''(t) to be 
r^'' = constant/wfc(t). These approximations amount to replacing the self-energies on the right-hand sides of ( 5.10| ) 
and (5.11) with constant equilibrium values. With the additional approximation that Tk <?C Wfe, we find 



2{5nk + Svk) 



peq 
^ k 



uJk^T^ 



1 + 2n'=i) 



(6.8) 



which reproduces the sum of ( |2l5| ) and ( |2^ ). 

It is not surprising that the above strategy agrees with linear response theory only when ujk- This indicates 

only that the methods available for resumming self-energies in the equilibrium and non-equilibrium theories differ 



13 



beyond the leading order in Vk/iOk- The fact that we can recover the linear response result only by ignoring the 
fluctuations in self-energies induced by those in the is, however, rather more significant, as we discover by attempting 
to solve the next-to-leading order equations without this extra approximation. To simplify matters, we continue to 
retain only the leading terms in Tk/(^k, in which case the order e equations are 



1 

2Ji 



-ujkdt (Imi/j,^) 



(6.9) 



dk'K{k, k')Snk' = dtnl 



where 



K{k,k') 



SSk 



Suk' 



'^k-r, eq 

l^k 



(6.10) 



(6.11) 



While the first of these gives Si'k explicitly, the second is an integral equation to be solved for Suk- It turns out 
that this equation has no solution. The reason is that the scattering processes described by Sk conserve both particle 
number and energy. One easily finds that this implies the two sum rules 



dkPK{k,k') 



dkk^ujkK{k,k') = 



(6.12) 



valid for all k' . The source terms on the right-hand side of ( 6.10| ) do not respect these sum rules, so the equation is 
not self-consistent and has, in principle, no solution. At higher orders of perturb ation theory, particle number is not 
conserved; only the energy sum rule remains, but that is sufficient to invalidate (|6.10| ). 

It is important to emphasize that the evolution equation (3.5) with e = 1 is perfectly sound: it is a Boltzmann 
equation with a source term, which presumably has a satisfactory solution for nk{t). What we have found is that this 

solution does not have a time-derivative expansion. That is, it cannot be expressed as nk{t) — nf\(l)) 



Nor, therefore, can the equation of motion (1.5). Our principal conclusion, then, is that the friction coefficient r]((j)) 
does not exist. As it stands, this conclusion rests on an approximate treatment of a particular model, the theory. 



It is likely, however, to be quite generic, as we discuss in section IX 



VII. NUMERICAL INVESTIGATION 



Although we have just reached the conclusion that the friction coefficient ?7((/)) does not exist, we have also seen that 
the linear-response result for ?7(0) can be recovered by ignoring fluctuations in the self-energies — an approximation 
that, at first sight, would seem not to be severe for a system reasonably close to equilibrium. Thus, although ri{(j)) 
formally does not exist, the local equation of motion (1.1) with ri{(f)) as given by linear response theory might be a 
reasonable approximation to the true equation of motion. We have obtained numerical results that may bear on this 
question by taking advantage of the following circumstance. A discretized and truncated version of (6.10) that one 
might attempt to solve numerically is 



KkM'Suk' = bk, 



(7.1) 



k'=0 



where bk stands for the source terms on the right-hand side. Because the kernel K k,k' n ow involves only values of k 
and k' up to the cutoff value A:max, it does not exactly obey the sum- rule constraints (6.12) and the truncated equation 
may have a solution. 

In fact, we find that it has a very well defined solution, as illustrated in figure ^, where the quantity plotted is 
Sx = {4:m^c/X(j)(j))Sn and we define the naturally-occurring coupling constant c — A^7r^/64(27r)'^. The example shown 
in figure ^ is for (3m^ = c = 1. The kernel Kk^k' decays rapidly for k ^ k' , but is not small when k ^ k' . Thus, 
when both k and k' are bounded by the cutoff k^n^, the sum rules may be verified for k' <C fcmax, but they fail 
for larger values of k' . Figure ^ illustrates this for the energy sum rule, with (3mip — 1 and a cutoff fcmax = 3Qni,^. 

The circles in this figure represent the integral J^^""^"^ dkk^ujkK{k,k'), while the squares show, for comparison, the 



quantity 0.1 J^^'"^'' dk k^ujk\K{k, k')\. The source bk in (p^.lj) becomes very small when k is greater than a few times 
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FIG. 2: The solution of (7.1) for /3m^ = 1 and c = 1 
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FIG. 3: Illustration of the energy sum rule (3.12); the quantities plotted are explained in the text. 



rritp, and figure 2 shows that the same is true of the solution Suk- In effect, we see that, regardless of the cutoff, only 
a "self-truncated" kern el, with k and k' restricted to values smaller than a few times m^, contributes significantly 
to the solution of (7T). The fact that this "self-truncated" kernel does not in itself respect the sum rules accounts 
for the existence of a well defined solution and, because of the self-truncation, we are able to verify that this solution 
converges to a cutoff-independent form as fcmax is increased. 



Formally, we can use this numerical solution in (4.11) to obtain an estimate for the friction coefficient 77. The result 
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FIG. 4: Friction coefficient calculated from the solution of ( [7.l| ) 



this for a range of coupling strengths is shown in figure ^, where the quantity plotted is u = rjm^p/ (1287r^^. 
shows the ratio rj/rji^^, where t^lr is the friction coefficient calculated in linear response theory from ( 2.1i ) 
and (2.21) in the limit Tk <^ cok- (The haphazard appearance of this figure results from the different dependencies 



of doin; 
Figure 



of rji and rju on temperature and coupling strength.) The formal quantity shown in figure |4| has negative values at 
weak coupling (where our perturbative methods are most likely to make sense) , and clearly cannot be interpreted as a 
genuine friction coefficient. This, of course, is consistent with our earlier argument that the friction coefficient is not 
well defined. The conclusi on is that non-equilibrium methods are needed to investigate the real effect of dissipation 
in the equation of motion (1.3), even for a system that may be quite close to equilibrium. One approximate m ethod 
of doing this is to integrate this equation of motion simultaneously with the evolution equations ( 4.16 ) and ( 4.17 ) 
and we plan to report on such calculations in future work. The large discrepancies apparent in figure |^ suggest that 
quantitatively significant deviations from the predictions of linear response theory may be expected. 



VIII. FRICTION ARISING FROM FERMIONS 



We comment briefly on the frictional effect of a Yukawa coupling to fermions. Non-equilibrium perturbation theory 
for fermions is discussed in Ref. [E7|, to which we refer the reader for the somewhat cumbersome details. With 



approximations analogous to those described above for scalar fields, we find for the relevant term in (1.3) 



l) + 2\k\Reiy, 



(8.1) 



where n\1'\t) and v^f^\t) are the fermionic analogues of the functions {t) and (t) 



the form 



{t), then the on-shell processes if 



If the fermions have a mass 
ipip are kinematically allowed. We then obtain kinetic equations of 



Otn):' = 



dtv];' = 



-r 



-2i 



Rei^ 



Ifcl 



OJ 



- 1) 



(8.2) 
(8.3) 



4'^)andr('^) 

now acquires a contribution from a fermion loop. Taking this into account, the kinetic equations (4.16) and (8.2) are 



where ai'^' and T]!^' are extracted as in ( B.lOj ) and (5.11) from the fermion self-energy. Naturally, the scalar self-energy 
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FIG. 5: Ratio of the formal friction coefficient of figure ^ to that calculated from linear response theory. 



a consistent pair of Boltzmann equations, in which the scattering integrals preserve the particle number 2N^ 
(we have derived 



2 ) and ( |8.3D only for zero chemical potential, in which case fermions and antifermions are equally 



abundant and means the total number of these particles) and the total energy of scalar and fermionic particles. 
These conservation laws (or, at higher orders, just energy conservation) again imply that neither the solution of the 
Boltzmann equations nor the equation of motion for <j) has a time-derivative expansion. 



IX. DISCUSSION 




The equation of motion ( |1.3| ) is inherently non-local in time; it represents a non-Markovian process in which 
evolution depends on the history 0(t') at all times prior to the time t of interest. When the state of the system is not 
too far from thermal equilibrium, it is tempting to suppose that a local equation of motion equation (1.1) might be 
approximately valid, the friction coefficient being estimated from equilibrium statistical mechanics. In this paper, we 
have examined the approximations needed to extract a local equation of motion from the non-local one, and concluded 
that this cannot in fact be done. Under suitable conditions (the principal requirement is the existence of a relaxation 
time short enough to ensure that co rrela tions decay rapidly compared with the characteristic ti me sc ale o n wh ich 
(j){t) changes) the expectation values in (1.3) can be approximated by local expressions of the form 
which the auxiliary functions nk{t) and viJt) themselves obey local kinetic equations, such as (4.16), 
( ^.31 ). However, this set of local evolution equations can be reduced further to a single local equation for (t>{t) only if 
the kinetic equations admit a solution in the form of a time-derivative expansion — and we find that they do not. The 
obstruction arises from fluctuations in self-energies, of which the equilibrium theory takes no account. 

Although our explicit computations focussed on the simplest example of a single, self-coupled scalar field, we 
have indicated is sections | and VIII that the situation is quite generic. The above conclusion emerges from an 
approximate treatment of non-equilibrium dynamics, which can hardly be regarded as a rigorous proof. It would 
seem that some approximation more or less equivalent to the local ansatze (3.27) and ( 3.2S| ) for self-energies is an 
inevitable step towards the derivation of a local equation of motion; without some such approximation, the expectation 
values in (1.3) remain non-Markovian and, a fortiori, cannot be represented by a local friction term. A subsidiary 
approximation made in section ^ was to set the quasiparticle width Tk to zero for the purposes of estimating self- 
energies. Although this approximation greatly simplified our analytical analysis, it can and should be avoided in a 
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comprehensive numerical study. The effect of this approximation is to restrict the scattering processes in Boltz mann 
equations to on-shell processes. Now, the inclusion of off-shell processes might well invalidate the sum rules ( |6.12 ) from 



which we concluded that the friction coefficient 77 does not exist. Formally, then, by including off-shell processes, we 
might after all be able to obtain a time-derivative expansion of the equation of motion. However, the friction coefficient 



implied by this expansion would be quantitatively similar to that obtained in section VIl. At weak coupling, it is 
negative (and thus physically unacceptable) and quite different from the one yielded by linear response theory. At 
strong coupling, the perturbative methods employed both here and in linear response theory are quantitatively, and 
perhaps also qualitatively, unreliable. Our practical conclusion, then, is that even for a system quite close to thermal 



equilibrium the local equation of motion (1.1) does not furnish an reliable account of dissipation, whether or not 
off-shell processes serve to recover a formal time-derivative expansion. A thorough numerical investigation of the non- 
equilibrium evolution is therefore essential. At weak coupling, a numerical implementation of the evolution equations 
developed here is quite widely applicable, and may well be quantitatively adequate, though recently developed methods 
based on the 2PI-1/N formalism |2^, are probably more powerful in situations where they can be applied. 

Finally, we observe that the analysis given here applies to quantum field theory in an expanding spacetime with 
only minor modifications. In a spatially homogeneous Robertson- Walker spacetime, the field redefinitions $ a^^$ 
and ■0 a~^/^V': where a{t) is the scale factor, serve to cast the theories we deal with in the form of a Minkowski- 
space theory with time-dependent masses, provided that t is taken as the conformal time coordinate. In the case of 
a spinor field or a conformally coupled scalar field, these masses are given simply by m(t) — a(t)m. Consequently, 



the formalism we have constructed changes only insofar as the masses in (1.4) depend on t through both (j){t) and 



a{t). This additional time dependence modifies evol ution equations in a way that may be cosmologically important. 



However, its net effect on, say, equations (6.9) and (6.10) is just that the right-hand sides of these equation contain 
terms proportional to d in addition to those proportional to 0. These induce additional contributions to 6nk and 6vk, 
but do not affect our conclusions concerning the terms proportional to 0. 
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